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Abstract 

A perturbation method deals with dipolar coupling spins in solids is presented. As example of 
application the method, the multile-quantum coherence dynamics in clusters of a linear chain of 
four nuclear spins and a ring of six spins coupled by dipole-dipole interaction are considered. The 
calculated OQ- and 2Q intensities in a linear chain of four nuclear spins and 6Q intensity in a ring 
of six spins vs the duration of the preparation period agree well with the exact solutions (for linear 
chain of four nuclear spins) and simulation data (for linear chain of four nuclear spins and a ring 
of six spin). 



Introduction 

The dipole-dipole interaction (DDI) is widespread in nature and determines behavior of 
many physical systems. In solids state nuclear magnetic resonance (NMR), it is respon- 
sible for many specific phenomena and plays an important role in nuclear spin dynamics 
^|. However in solids the evolution of a spin system under the DDI involves many spins 
and behavior of such system cannot be analyzed analytically. Even the numerical analysis 
becomes difficult because the number of states = 2" is growing exponentially with the 
increasing of number of spins n. Therefore existing theories are based on phenomenological 
models and only macroscopic characteristics such as spin-spin relaxation times, the second 
and the fourth moments of resonance lines are taken into account [Ij. These difficulties are 
very clearly displayed in multiple-quantum (MQ) spin dynamics. The MQ phenomena in- 
volve various multiple-spin transitions between the Zeeman energy levels and MQ coherence 
is formed at time t > tu^^ , where lu^ is the characteristic frequency of DDI ^Q|. Hence, 
ujdt > 1 is not a small parameter, that did not allow one to use perturbation methods to 
study MQ dynamics. Indeed only simple exactly solvable models of a spin system such as 
iwo and three dipolar coupling spins 1/2 Q| or one-dimensional linear chains of spins 
isl ] were analyzed theoretically. The last achievement in this direction is the model with 



15l |. Note that the simplified calcula- 



identical DDI coupling constants for all spin pairs 
ions essential for the case of identical DDI coupling constant have been already mentioned 
irj l . This model describes only zero- and second-order coherences. But in real solids, DDI is 
characterized by the different coupling constants. In order to obtain important information 
on molecular structure and spin dynamics of these systems, the evolution of large number 
modes of MQ coherence has to be analyzed. 

Another important purpose of the study of MQ coherence is the field of quantum infor- 
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mation processing p, lD, |8[. It was experimentally demonstrated the possibility of creating 
pseudopure spin states in clusters of coupled spins by MQ method |3] and that dynamics of 
;he quantum entanglement is uniquely determined by the time evolution of MQ coherences 

Our objective is to adapt a perturbation approach to the problem of dipolar coupling 
spin dynamics in solids and to obtain the description of large number modes of MQ co- 
herence. The Hamiltonian of the spin system with different DDI constants is divided into 
several parts, each of them is characterized by its constant. Our main idea is to take into 
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account in MQ dynamics the influence of tlie different items of the Hamiltonian with dif- 
ferent degree of accuracy. The items corresponding to smaller DDI coupling constants are 
considered as a perturbation while the item with the largest constant (describing the inter- 
action between nearest neighbours) is taken into account exactly. The proposed approach 
will be a power method to describe wide range of physical problems dealing with dynamics 
of dipolar coupling spins in solids. On the one hand, this approach uses the advantages of 



On the other hand, it simplifies calculations by using a 



exactly solvable models 
perturbation technique. 
Theory 

Let us consider a system of nuclear spins 1/2 placed in an external magnetic field i^oll ^ 
axis and subjected to a suitable designed pulse sequence, thus permitting the observation of 
MQ coherences The effect of the sequence of irradiating pulses on the spin system can 
be represented by a unitary transformation propagator [l^ 

f/(t)=e-*^, (1) 

where Ti is the effective time-independent Hamiltonian. First assume that the effective 
Hamiltonian can be divided into two parts with only two different DDI constants a and (3: 
Ti = A + B , where the operators A and B can be presented in the following form: A = aA 
and B = [3B, a = \\A\\ and /3 = are norms of operators A and B, respectively (a > /3 
and [A,B] ^ 0). The operators A and B do not include any DDI couphng constants. We 
present expression (1) in the form 

e-**("^+'^^)=e-**^^aA(t), (2) 
where operator a {t) obeys the differential equation 

t^ = aA{t)aA{t), (3) 

with the initial condition 

CTA (0) = 1 (4) 

and A{t) = ^-^^^P^ j^^itpB ^ Assume that a > (3. Solving (3), we limit ourselves by the first 
power of -. Expanding the exponential term in the right hand of Eq. (2) in a series and 
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substituting in Eq. (3), we obtain 

i^^^=a{A-itf3[B,A])aA{t). (5) 

To solve Eq.(5) we will use the iterative method. We will search the solution of Eq.(5) as a 
series on parameter (at)"' : 

CTA(t)=J2^^A\t) (6) 
n=0 

Taking into account that cr^^ (t) — 1 the solution of Eq.(5) for n = 1 can be 
obtainedia^ (t) = {yA - \y'^e [S, A]) 

(t) = {-iat) A. (7) 

For n — 2, ones obtains 



<'{t)^'-^[A' + ^[B,A]j (8) 
Keeping only linear in ^ terms the following expression for operator cr^'* (t) can be obtained 



^ { tatr -^(^{m- 1) BA^-' - A'^'^-^BA^ j (9) 

Utilizing (9) to expand the expression of the unitary transformation propagator (1) we 
find 

[3 ^ {tatf^^ ^ (-1)^' 



U (t) = e-^*("^+^^) = 1 _ ^ y i!^^ y , ~^ ,, A^BA^-^ e— (10) 
Taking into account that 

y TTT— = 4 {B, (11) 
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ill 

j=0 

where the repeated commutators bracket {...} is defined by 



{B, A'} = B, {B, A^+'} = [{B, A^} , A] , (12) 

we obtained that expression of the unitary transformation propagator takes a form 

U (t) = e-**("^+/^^) = h - ^ ^ {B, j e-^"*^. (13) 



The summing over n in (13) results in 



= ^« / f^te'"*^Se-*'^*^. (14) 

n=0 J 



Substituting Eq. (14) into (13) we obtain 



latA 



(15) 



which is a well-known formula for expansion of an exponential operator in a perturbation 
series jiol. l2o|. Expanding the operator exponent e~*°*^ of the left-hand of Eq. (13) gives 

^-uiaA+m = y (z!^ (a'^-^Y {B, A-} A^ I (16) 

^ k\ \ {n + l)\ ^ ^ I ^ ^ 

k=0 \ n=0 ^ ^ 

After rearranging of the series (16) we obtain 



^-it{aA+f3B) 



a ^ ^ n\m\ (n + m + 1) 

n=0 m=0 ^ ^ 



g-*atA (17) 



In the same way we obtain the expansion up to second order in the ratio 



^ oo oo / 1 \m / • ,\n+m+l 
ji=0 m=0 ^ ^ 

a/ '^^n!m!A;!/!(A; + / + l)(A; + / + m + n + 2)^^ 

Expression (18) can be generalized to a case when the exponential operator contains 
arbitrary number of the non-commutative operators and can be extended to include various 
powers of the operators. For example, if the effective Hamiltonian operator includes three 
non-commutative operators 'H = A^E^C = aA + [3B + 5C, where 5 < /? < a, we can 
obtain the expansion in series for a three-spin cluster replacing B hj B + |C in Eq. (18) . 

Expressions (17) and (18) represent independent expansions in two parameters £ < 1 
and X = at. Eqs. (17) and (18) appear to be complicated but in fact it is quite simple to 
use, as the following examples will illustrate. 

Results and discussion 

Let us consider a cluster which constitutes a linear chain of four nuclear spins 1/2 coupled 
by DDI. The MQ dynamics in the rotating frame is described by the time-independent 
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average Hamiltonian: 

^ = E ^^-^ = E {It It + ijik) (19) 

j<k j<k 

and it and /~ are the raising and lowering operators for spin j. The dipolar coupling 
constant, djk, for any pair of nuclei j and k, is given by 

d,-, = ^(l-3cos'^,fe), (20) 

where 7 is the gyromagnetic ratio of the nuclei, rjk is the internuclear spacing, and 9jk is the 
angle that the vector fjk makes with the external magnetic field. In the high-temperature 
approximation the density matrix at the end of the preparation period is given by 

= e-^^V(0)e'^* (21) 

where p (0) is the initial density matrix in the high-temperature approximation 

P(0) = E^I' (22) 

Ij is the projection of the angular momentum operator of spin j on the direction of the 
external field. The average Hamiltonian (19) can be divided into the three parts according 
to the number of different coupling constants di2 > dis > du : 



H — H12 + Ti-ia + "^14, (23) 

where di2 is the coupling constant of the nearest neighbors, ^13 is the coupling constant of 
the next nearest neighbors, and du is the coupling constant of the next-next-nearest neigh- 
bors. The experimentally observed values are the intensities of multiple-quantum coherences 

-^^Q (^) ^ j^j.p2 (0) E pig (^) n = m,p- m,g. (24) 

where m^p and m^g are the eigenvalues of the initial density matrix (22). Using expansion 
(18) up to second order in ^ = ^ and f = ^ and keeping terms up to eighth order 
in X — d\2t , one can determine the normalized 0-quantum (Jqq) and 2-quantum {J2Q) 
intensities: 
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(25) 



(26) 



In the case where the parameter x is not small we have to take exactly into account of 
the terms containing x . By summation of series in (12) over n, m, I, and k we obtain 
exact analytical expressions which give the time dependence of intensities of - quantum 
coherence 



Jqq — 1 ~ 2 



2Q 



(27) 



and of 2-quantum coherence 

Jon = — T + T cos(a;\/5) cosfa;) ^ sm(xV5) cosfa;) ( -\ + 3 ( — 

^ 4 4' 50 \a J \a 

2 J cos(a;\/5) cos(a;^ 



—X 



a J 4 V a 



in which the interaction between the nearest neighbours is taken into account exactly. Now 
let us compare the time dependence of the MQ intensities given by (25) - (28) with those 
obtained using an exact solution of the evolution equation and using computer simulation. 
The exact solution gives the following expressions for Jqq 



OQ 



1-2J: 
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(29) 



and for J2Q 



J, 



2Q 
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(30) 
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where 



K±={-]±1 (31) 



a 




P= ^ +2 r +5 + 4^ (32) 



P± = R+± R. (33) 




-2(^)+5±8^ + 4^^J^j (34) 

Figs. 1 and 2 show the evolution of the normahzed OQ and 2Q coherences in the hnear 
chain of four nuclear spins 1/2 coupled by DDI, where £ = | ? f = ^ and at t = the spin 
system is in thermal equilibrium determinates by matrix (16). These figures present also the 
results of computer simulation performed using the MATLAB package. One can see that 
expressions (25) and (26) gives a good agreement with the exact analytical and numerical 
results for X < 1 , while expressions (27) and (28) agree up to x = 5. Note that our results 



coincide also with the dependences obtained in paper 

As a second example let us consider the MQ dynamics in a ring of six spins coupled by 
dipole-dipole interaction. The ratios of the dipole-dipole coupling constants are given by 
a ~ 373 ^^"^ a ~ I' "^here a is the coupling constant of the nearest neighbors, /3 is the 
coupling constant of the next-nearest neighbors, 6 is the coupling constant of the next-next- 
nearest neighbors. Calculation using expansion (18) up to second order in ratios of £ and 
^ and eighth order in x = at gives zero intensity for 6Q coherence. The non-zero intensity 
can be obtained by taking into account exactly of terms related to the interaction between 
the nearest neighbors . Using Eq. (18) we obtain the analytical expression of the intensities 
of 6Q coherence calculated up to second order in ratio of £ and ^ 

JdQ = - ,r.^L. [84a; (2 cos x + cos 2x) + 88 sin x - 5 (32 sin 2x + sin 4x)? (35) 
497664 

It is interesting that Jqq coherence up to second order does not depend on the interaction 
between the next-nearest neighbors. The validity of the time dependence of Jqq according to 
Eq. (35) was tested by comparison with numerical results. Fig. 3 shows that the prediction 
of Eq. (35) gives a reasonable agreement with the simulation data up to x < 7 , which 
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coincides with the estimation using the dipolar couphng constant: x < f = 8. The results 

nn 

of our simulation, shown in Fig. 3, coincide with the dependences obtained in Refs. |13l.ll8|. 
Conclusion 

In conclusion, the perturbation method was developed which is based on the expansion 
of the operator exponent in a perturbation series. It allows us to apply the perturbation 
approach to the description of the MQ spin dynamics in solids. The exact and perturbation 
analytical expressions were obtained which describe OQ and 2Q dynamics in a cluster in the 
form of a linear chain of four nuclear spins coupled by DDI in solids. The calculated OQ- and 
20 intensities versus the duration of the preparation period agree with the simulation data 
j2l|. The obtained analytical expression for 6Q dynamics in a ring of six spins coupled by 
DDI is in a good agreement with the results of numerical simulation. The developed method 
can be extended to include various power of the operators and applied to the description 
of a wide range of physical problems that deal with dynamics of dipolar coupling spins in 
solids. 
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Captions for figures. 

Figure 1 

Fig. 1. Time dependences (in units of ^ ) of the normalized intensities of OQ coherence 
in a linear chain of four nuclear spins 1/2 coupled by DDI. Solid line is the exact solution 
(Eq.(29)), dotted line is the calculation using Eg. (25) , da shed hne is the calculation using 
Eq.(27), open circles present the simulation results of [2l|, and down solid triangles are our 
numerical simulation. 

Figure 2 

Fig. 2. Time dependences (in units of ^ ) of the normalized intensities of 2Q coherence 
in a linear chain of four nuclear spins 1/2 coupled by DDI. Solid line is the exact solution 
(Eq.(30), dotted line is the calculation using Eq.(28), dashed line is the calculation using 
Eq.(26), open circles present the simulation results of and down solid triangles are our 
numerical simulation. 

Figure 3 

Fig.3 

Time dependences of the intensities of 6Q coherence in a ring of six spins. The dashed 
line is the prediction of Eg. (35) up to second order in ratio of and the solid line is the 
simulation data from Q • 
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